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Summary 


A theory of non-linear crankshaft vibrations is developed, 
fully considering the effects of impulse distortion and mass 


variation. A case of "secondary resonance" in a 10-cylinder 
engine is described and calculated. Representing the non- 
linearities as "fictive forces" acting on a linear system a great 


simplification of detail calculations is obtained. 
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Notations used in the paper 


a crank angle, measured from the "main position'' where the 
total moment of inertia is minimum. 

4 2 

28" 
B. = coefficient of oscillating mass pr. cylinder, 

c 
Bac = BA(- x, 1, A) coefficient of mass variation pr. cyl. 
of order p = (1, 2, 3). 
A) 


= coefficient of total mass variation in 


multi-cylinder engines, 


7 short for 
A ) 
= coefficient of coupled mass variation (of 
2m 4, 2. order) from vibration form A to B. 
Y phase angle of impulse at main position. 
- - - vibration at - 


angle between cranks. 


4 relative vibration amplitudes in a vibration form. 
= dimensionless coefficient of total damping. 
wow 
e 
= connecting rod relation. 
x vibration amplitude, radian, 


a, b, c maximum values of vibration amplitudes at the point 
where A= 1. 


a,ora = reduced or relative impulse of order p, 
P w.ZmA° also called "static amplitude", radian. 
2(K_ 4”) 
ay. reduced impulse for subharmonic excitation. 


wo 


a same, with oscillating-mass impulse included in x: 


a' _, b' etc, fictive force or artificial impulse of order indi- 
cated and produced by vibration denoted by the letter. 
Measured relatively as a, and a. 


F == magnification factor. 


8 
A) 
8 
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J total, equivalent inertia moment, variable, 
=m a® - Mean value. 
K, cylinder impulse of order p pr. revolution. 
q number of vibrations pr. revolution. 
w, angular velocity of crankshaft. 
w phase velocity of vibration. 
w - - - free vibrations, calculated for constant 
masses, 
w 9% 
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4. Introduction 


In certain types of 10-cylinder two-stroke engines torsional 
vibrations of 8. order have been observed, which could not be ex- 
plained by the common calculation methods, as the geometric 
resultant of the 10 impulses of 8. order for the adopted crank 
position was very small and should not be able to exite vibrations 
of a size more than about one fifth of the actually measured. 

Luckily, these cases were not of seriuus consequences as the 
vibrations were found somewhat above normal revolutions of the 
engine and would therefore not occur during normal service, but as 
such an unexpected phenomenon might, in other cases, arise at a 
more inconvenient range, a full explanation of it was wanted in order 
to avoid it. 

Thorough investigations showed, that there could be no consi- 
derable action of impulses from the propeller, coupled vibrations 
from sources outside the engine, etc. Furthermore, the stress curve 
about the resonance point was of a shape similar to ordinary re- 
sonance curves as if the impulse of 8. order were about 5 times 
bigger than it really was, and this artificial, extra impulse must be 
sought within the vibrating crankshaft system itself. 

Hereafter it was clear, that the vibrations must be a kind of 
"non-linear" vibrations which cannot be calculated by the common, 
"linear" vibration theory. Some knowledge of non-linear vibrations is 
thus necessary for those responsible for the calculation of torsional 
vibrations in a great variety of big and costly engines. Even if non- 
linear vibrations can be expected only in a few and special cases, 
the study of them is in any case interesting and will give an in- 
creased feeling of safety also in the normal calculations. 

The appellation "linear" is here used for systems determined 
by linear differential equations with constant coefficients, and ''non- 
linear" for systems with variable coefficients, whether the varia- 


tions are time functions or functions of the vibration amplitude x. 


In the mathematical terminology only the latter systems are truly 


| 
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"non-linear" (namely non-linear in x), but in vibration theory there 

is not much difference between the two kinds of variation because x 
is in itself a periodic time function when there are steady vibrations, 
and it is therefore more important to have handy names to distinguish 
between systems with constant coefficients and with variable coeffici- 
ents. Yet, when we have to distinguish between the two kinds of va- 
riation (the time-depending and the vibration-depending) we shall here 
call the former systems "freely variable'' and the latter ‘true non- 
linear systems", while the cases with both time- and vibration-de- 
pending variations (to which the crankshaft system belongs) will some- 
times be called "guided variable systems''. These names have a dis- 
tinct physical meaning, and the "freedom" or "guidance" of the va- 
riations will always appear directly from the constructive features of 
the system. 

The problems concerning a ''second order theory" of torsional 
vibrations were some 30 years ago thoroughly treated by Trefftz, Kluge, 
Grammel and others), and the subject may have been considered as 
exhausted with these highly theoretical papers. But while their funda- 
mentals, chiefly based on application of Lagrange's equations, are of 
unchanging validity, their subsequent treatment of technical details and 
approximations seem to-day unsatisfactory and insufficient from the 
engineer's point of view. On the other hand, there is to-day an exten- 
sive litterature on non-linear vibrations, but neither is this sufficient 
for the use of the mechanical engineer meeting new and unexpected 
problems in his practice. Especially there is very little to be found 
about non-linear vibrations in multi-mass systems, although dangerous 
vibrations in practice always concern multi-mass systems. It has 
so far been left to the practical man to find out how he shall best 
substitute his real, multi-mass system with an equivalent one-mass 
system which the mathematician will calculate for him with any 
wanted accuracy and subtlety. 

The investigation of the above mentioned unexpected vibrations 
resulted in a rather complete theory of non-linear crankshaft vibra- 
tions and easy and yet in first approximation completely correct cal- 


culation methods as described in the following. Besides, it was 


1) see references at end of paper. 
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necessary to develop certain basic principles for a general theory of 
non-linear vibrations in multi-mass systems. These matters are 
mostly found in sections 4 and 8-10. 

The parts printed with closer line distance contain (besides de- 
tail calculations for the control of energy balances etc.) the theory 
of subharmonic vibrations, while the main text goes straightly for the 
theory of ''secondary resonance'' which seems to be the more im- 
portant matter, in any case for crankshaft vibrations. A general theory 


of non-linear vibrations must comprise both phenomena but also clearly 


distinguish between them. 
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2. 


Linear vibrations in one-mass and multi-mass systems 


It will be appropriate first to recapitulate the calculation of 
steady vibrations in linear systems with constant coefficients, and 
to introduce a special terminology which is very convenient, not to 
say necessary, for treating the intricate problems of non-linearity. 


In a one-mass fixed system fig. 1, consisting of a swing 


Vibration form 
| stittness | ax 
Fig. 1. 


Torsional vibrations in a linear one-mass system, 


mass with inertia moment m kgcmeec” and a shaft with stiffness 
s kgcm/rad, acted upon by an impulse K sin(wt + Y) kgcm and a 
damping ro kgcm, the angular vibration amplitude x is deter- 
mined by the differential equation 


m SF +r = Ksin(ut + ¥) (1) 
t 


Dividing with mun = s we get: 


+ +x= a, Sin(wt + Y) (2) 
w dt e 
e 
where a, ="; means a reduced or relative impulse of dimension 


vibration amplitude (here dimensionless), also to be called "static 
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amplitude" as it is the deflection which the system would get for a 


static force K, and p = — is a dimensionless damping coefficient. 
e 
Introducing §— = , the relative frequency of the impulse in 


relation to the frequency of free vibrations, we get 


2 


4... dx 
+ pet + x = a_sin(wt + Y) 3 
d(wt (3) 


} The continuos solution (disregarding tuning-in) is 


x = asin(wt +9) where a= a.:F 
4 7 
the magnification factor F= 
2\2 2,¢ 
and the phase difference is given by , 
-@) = (4) 
or sin(y -@) = pEF 


It is very instructive to draw a vector diagram, fig. 2, for the 
forces expressed in equation (3) by their relative values. They are 
all considered as outer forces acting on the swing mass like the 
impulse, and the gotating vectors are shown in their position att= 0. 
From this diagram we easily get formulae (4). In the resonance 
case —& = 1 we have F =< and Y -9= i The impulse is leading 
the vibration by 4 pore Disregarding damping outside of reso- 
i-€ 
—& > 1. The impulse is here acting like an elastic force decreasing 


nance we have F = 


with Y - p= 0 for — < 1 and = x for 


or increasing the real elastic force and thereby the frequency. 

This is of course well-known, elementary vibration theory, 
but expressed in the simplest possible way, and a complete fa- 
miliarity with the terminology is necessary for the following. 
Furthermore, it is here of some special interest to note, that under 


steady conditions the vibration phase 9 will automatically adjust it- 


self to the value given in (4), which is easily seen from the vector 


| 

4c 

| 


diagram. If the vibration should 
lack behind, the resulting elastic 
force will be bigger, as a bigger 
component of a. will act as 
elastic force. Hence the frequency 
of the vibration will increase, and 
the vibration will advance until 9 
has obtained the correct value. 
Oppositely, if the vibration by 
some disturbance momentarily 


should be too far advanced in re- 


lation to the impulse, the fre- 
quency will decrease and the vi- 

Fig. 2. Vector diagram bration fall back to its right po- 

for steady vibrations. sition. Thus the phase difference 

given by (4) is a stable equili- 
brium. 

In a multi-mass, free system fig. 3, consisting of several 
swing masses on a freely rotating shaft, acted upon by a number 
of impulses of different size, frequencies and phase angles, the 
steady vibrations can always be analyzed in their components of the 
different free vibration forms of the system, of which one is shown 
in the figure. This is by far not so well-known as the elementary 
one-mass theory, but nearly indispensable for the correct calcula- 
tion of any complicated vibration system. 

Each of the said vibration components can be calculated 
exactly after the same formulae as for the one-mass system, when 


a. and p are specially calculated for each separate vibration form 
by means of the expressions 


2 
2 mA w 2m & 


All the masses are namely moving proportionally (in that par- 
ticular vibration component) with relative amplitudes A, and their 
total inertia is therefore the same as for a single mass mA 
placed at the point 1 where A is taken as unity. 2(KA) represents 
the geometric sum of the impulse work vectors (of the same fre- 


8 
| Inertia force 
2 
: af 
fa 
er force | \ w 
a 
| a impulse 
5 


Fig. 3. Linear multi-mass system, 


quency) and means a resulting impulse which, working at the point 1, 
will produce the same impulse work (on that particular vibration 
component) as all the real impulses of that frequency. Similarly, 
w2ra represents the numerical sum of the damping work vectors 
in the resonance case from which the total relative damping coeffi- 
cient p is derived by the same reduction as used for the impulses. 
Then the magnification factor and the phase difference between the 
resulting impulse and the vibration are, for that vibration component, 
calculated exactly after the formulae (4) set up for a one-mass 
system. 

The given impulses should in the same way be placed on each 
of the other possible free vibration forms of the system, and the 
calculated vibrations should all be added (geometrically) giving the 
real, total vibration, while the calculated vibration components in 


the separate vibration forms were only abstractions. This may be 


called the principle of formal superposition, generally not recognized 


in elementary handbooks although it is easy enough. 

The same set of given impulses is thereby used as many times 
as there are free vibration forms worth while to calculate, but after 
adding all the calculated vibrations we have in reality only used the 
full value of the impulses once. This follows from the “orthogonality” 
of the free vibration forms, i.e. the rule, that for two different 
vibration forms A and B in a linear system (and also in a non-linear 


system, "linearized'' as explained in section 8) we have always 
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=m 4,4 0 


b = 
Only for a = b the product sum has a value, namely 2m a? . 
When forming the 2(K4,) for a vibration form A we have in 


reality only used impulse components (on all masses) 


2(KA,) 
K, = mA, 

as these will give a total impulse work equal to the work of an 
impulse 2=(K4,) while all other impulse components of the form 
C- mA, would give an impulse work 0 on the vibration compo- 
nent A. But such a set of ''resonance impulses" K, proportional 
to m4, for each mass will always give vibrations of the vibration 
form A independently of the frequency. On the other hand, adding 
all the so defined resonance impulses working at an arbitrary point 
in the system, the resultant must just be the real impulse (possibly 
no impulse) working at that point, because each resonance impulse 
is fully used in its own vibration form, so that the total vibration 
will require the resultant of the resonance impulses, while the given 
real impulses are, of course, in a unique way producing the total 
vibration, and so the two sets of impulses must be identical. 

When the impulses are at or near resonance with one parti- 
cular vibration form, it is usually sufficient to consider this vibra- 
tion component only, calculated by formulae (4) and (5), but it is 
of importance for the following (and in many other cases) to know, 
that the method is no accidental approximation but in fact com- 
pletely correct, if we add the vibration components calculated in 
the same way for all the other possible free vibration forms of the 
system (including the free parallel movement, A = 1 for all masses, 


= 0, , disregarding damping). 
2m 


= 


3. Non-linear crankshaft vibrations in general. 


Subharmonig vibrations due to impulse distortion 


In the differential equation (1) for a one-mass system, which 
may now also represent a single vibration component in a multi- 
mass system, there are 4 coefficients, m, r, 8, and K, which 
might be variable either with time or with the vibration amplitude, 
each of which possibilities would cause different kinds of non-linearity. 
But in crankshaft vibrations a variation of r, 8s, or K does not come 
into consideration in this connection. r is always so small that a 
variation of it cannot have any observable influence on the vibrations; 
s is constant as the vibration stress is always well within the pro- 
portionality limit of the shaft material; and K (meaning the maximum 
value of the impulse) may be variable with the load on the engine but 
not with time or vibration amplitude at steady load. 

Consequently, we should normally have to do only with variation 
of the swing mass m. Such a variation is always present and of 
considerable size pr. cylinder due to the variable contribution from 
the oscillating masses of piston and connecting rod. In practical 
calculations this variation is usually disregarded and the oscillating 
masses incalculated as a mean value, but this is in reality a very 
rough approximation although sufficient for most practical purposes. 

It is therefore our main purpose here to investigate all consequenses 
of the mass variation which might lead to important deviations from 
the common linear theory if only in special cases, 

If we now would ask whether this system with variable mass 
mathematically speaking is a non-linear system or a linear system 
with variable coefficients (a freely variable system), the answer 
must be, that freely variable is it only with some approximation, 
because the inertia moment is depending on the crank angle a, and 
this is not a pure time function. If the engine is running at a constant 
number of revolutions 4 the angle of rotation is wt at the node 
point, while the crank angle is a = wit + x. The crankshaft system 


is therefore,also in mathematical language, a non-linear system 


although the variation is mainly a time function but it is to some 
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extent also depending on x, and in our present investigation we 
cannot allow ourselves to disregard such a dependence without 
having calculated its importance. We should rather expect this 
dependence to have some real importance because it represents a 
physical reality, namely that the variation is ''guided" by the vi- 
brating shaft and not freely variable. 

Looking now once more at equation (1) to see whether there 
are other things depending on the parameter a, we find another 
source of non-linearity, namely the impulse itself. It is calculated 
by analysis of the tangentialpressure diagram as function of a, and 
it is solely depending on a as the cylinder pressure is given by 
the position of the piston. The cylinder impulse K, of order p 
periods pr. revolution is therefore in reality: 


K ,sin(pa +Y)= K +x) +Y) K 
(6) 
+ pxK ,cos(pw,t 


and there is an “artificial impulse’ which might be called non-linear, 
as its size is varying with the vibration amplitude, but if we place 
it at the left side of the differential equation, it is only a "linear" 
elasticity with variable coefficient (Matthieu's equation) - which shows 
better than anything else, that a variable linearity might as well be 
called a non-linearity. 

We have, accordingly, in the crankshaft system two sources of 
non-linearity, namely impulse distortion and mass variation. The 
latter will be our main subject, but we shall first look a little at 
the former. 

Although the artificial impulse (6) is always very small, as 
px << i, it has a general effect of some importance, namely its 
ability in the resonance case to pick up energy from the pure ro- 
tation of the shaft and transform it into vibration energy. Conside- 
ring first a one-cylinder engine it is seen, that when the impulse 
K sin(pa +Y) is in resonance with the frequency of the crankshaft, 
the very vibration which it excites itself, x = - acos(put +Y), will 
produce an artificial impulse: 


a 
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2 4 
paK cos = paK,, cos(2pw.t + 2y) (7) 
The new impulse of order 2p is usually without importance, 
only producing small forced vibrations, but the constant "impulse", 
- 5paK,, means a braking momentum which, when surmounted by 


the power of the engine, will pick up a work 
pr. rev. (8) 


equal to the work done by the main impulse Ksin(pw,t + Y) on the 
vibration - a cos(pw.t + Y) during p vibrations, and of course also 
equal to the damping work which is 


® praw mes = epaa_* mw xpaK 
e e 8 e P 

The cylinder impulses can excite the vibrations but cannot 
themselves supply vibration energy, as they are merely acting like 
elastic springs. But through the impulse distortion a braking of the 
shaft and an energy supply to the vibrations are effected at the same 
time, because the impulses considered as springs are not worked 
from outside, but from the rotating shaft itself. The impulses are 
"guided" by the vibrating shaft. 

Considering next a multi-cylinder engine at resonance in vibra- 
tion form A with the cylinder impulses of order p, the vibration 
synchronous for all cylinders may be a° sin(pwt +) at cyl. 1 and 
x=a Asin(pwt +) at the other cylinders, while the impulse is 
K sin(pa +Y) at cyl. 1 and K sin(p(at @) +¥) at another cylinder 
with crank angle © in relation to cyl. 1. It is then seen, that while 
the impulse work vectors for the main impulses are of a size pro- 
portional to x4 and situated at relative angles p®, the artificial 
impulses are already in themselves of a size proportional to K, 4 
and situated at the same relative angles. But the constant part of 
the artificial impulses (7) are braking momentums and should be 


added directly, consequently by means of the same vector diagram 


as used to calculate a_ = 2% . Therefore, as the total impulse 


3 
ag 


i44.~ 


acting on the system considered as one- mass system is 


the correct energy supply. 


a: =m a the total braking will be -$ paa, ° =m giving 


It is seen that any impulse of order p will (in multi- 
cylinder as well as one-cylinder engines) under the influence 
of a big resonance vibration of order q give artificial impulses 
of order p-q and ptq, but normally they will not themselves be 
in resonance with the vibration and will therefore be very 
harmless. The only dangerous case is 


P-q=q p = 2q 


in which case we may get a subharmonic vibration, a resonance 
vibration excited by an impulse of double the frequency of the 
vibrations. 

We shall first consider this case in a one-cylinder engine. 
It is convenient here to use the differential equation in the 
form (2). In stead of a, we shall, however, write ay meaning 


exactly the same for a one-cylinder engine, namely » but 


m 
for a multi-cylinder engine it will be found different from ag, as 
we shall see. For a one-cylinder engine with a relative impulse 
a, or ag of order p = 2q and in resonance with q vibrations pr. 
rev. we have then: 


2 


d“x dx = a_sin(2q(w_t + x) + Y) 
+x = a sin(pa + Y) u (2q(, ) 
e 


a sin(2wt + Y)* cos2qx + a cos(2w.t + ¥)* sin2qx 


(9) 


a sin(2w t + Y)° (4-2q7x") +a 2qx: cos(2w,t + Y) 


2 
= a sin(2w,t + Y) + 2qa x: cos(2w,t + Y)- 2q’a.x sin(2w_t + Y) 


this time using 3 members in the series for the sinus of a big 
plus a small angle, so that we have both a primary and a 
secondary artificial impulse, 

Introducing x = asin(wet +) in the primary artificial im- 
pulse (2. member in the last line of eq.(9)) we get from it a 
part of order q, that means an impulse of angular velocity 
G% = We, in resonance with the vibration. This impulse will be: 


a = - sin(w,t +Y-@) (10) 
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It is essential that the phase of the artificial impulse is 
depending on the vibration phase 9 but with negative sign. 
If it had been positive sign, the phase difference between im- 
pulse and vibration would have been unchanging, but as it is, 
the vibration is free to adjust itself to any suitable phase 
difference Y - 9 - @ = Y - 2@ between impulse and vibration, 
just as in ordinary resonance case explained in connection 
with fig. 2. In the present case eq. (9) concerns strict re- 


sonance, £ = 1, with — meaning —? , and the said phase 
difference will tend (as the impulse (10) is negative) towards 


in which case the impulse is directly opposed to the damping 
force, giving the condition for a steady vibration: 


pa=qaia P= qa, (14) 


While in usual resonance the condition pa = a, gives a 
determination of the amplitude a, the corresponding condition 
(11) in the subharmonic case leaves the amplitude underdeter- 
mined and yields in stead of it a value for the damping which 
will allow steady vibrations of arbitrary size. If the damping 
is bigger than given by (11), the vibration cannot come up, 
and if the damping is smaller, the vibration will immediately 
arise and increase to the utmost limits or breaking-down of 
the system. 

. lf there were no damping, p = 0, the condition for a 
steady solution would be that the impulse (10) is either in 
phase with the vibration or directly opposed to it, Y- 29=0 
or ®, in which case (10) would be an elastic force, either 
increasing or decreasing the frequency, so that we have in a 
way a "'splitting-up"' of the proper frequency in 2 resonance 
points (numbers of revs.) where steady vibrations can exist 
of arbitrary size, but between them there will be a range 
with increasing vibrations (in the mathematical language 
called "instability" - but nobody would expect "stability" at 
resonance when the impulse is bigger than the damping). 

Writing (9) in the form (3) we find the limits of the 
dangerous range given, analogous with (11), by 


2 4 
daa, 


In the general case with a damping p < qa, we get 


-29"2-s set 


a result which is also — by putting the ordinary magnifi- 
1 

cation factor (4) equal to ——— ee" " Gec° F is namely the mag 

nification factor for steady vibrations and in the present case 

valid only, when there are steady vibrations, 

We have calculated the secondary artificial impulse in (9) 
in order to show, that as the vibration is here excited by the 
primary artificial impulse, it is the secondary that will pro- 
duce the braking momentum supplying the energy. Introducing 


xfa sin(w.t +3 + 3) in the last member in (9) we find a 


constant part 


q 
and a braking work pr. rev. 


equal to the work pr. q vibrations of the impulse qaKp and 
also equal to the damping work in case of the equiltbriam con- 
dition (41). 

Comparing with (8) it is seen that the impulse work in 
subharmonic exitation is much smaller due to the extra factor 
a << 1, but at the same time of a much more dangerous 
nature as it is increasing with 2. power of the amplitude like 
the damping work, and the vibrations are therefore not limited 
by damping, if they can arise at all. 

The differential equation (9) and consequently also the main 
results (11) and (42) will be correct also for multi-cylinder 
engines when the relative value a, of the forcing impulse (of 
order p = 2q) is calculated as 


=(K a ) 
ae (13) 


in which 2(Kpd” ) means the usual vector summeatipa for vectors 
= order p and relative angles p®@ but of size K,A® in stead of 

K,4 as used for ag. The primary artificial impulses from the 
different cylinders are namely of the form x* Kpsinp(a+ ®), and 
as x is synchronous for all cylinders, it cannot change the rela- 
tive phase angles, only reduce the size with the factor A. 
Consequently the impulse vectors are to K,A and the 
impulse work vectors proportional to KpA’, but at the same 
angles as the ordinary impulses, For : point-symmetric vibra- 
tion form =(KA*) is very different from 2(KA), and the order 
numbers giving the biggest impulse summation for subharmonic 
excitation may therefore be different from those giving the big- 
gest ordinary excitation. 


= 


In practice the possibility of such vibrations is small 
(otherwise we should have heard more about them!), because 
qa, will usually be smaller than the coefficient of natural 
damping inherent in the crank mechanism itself, which is 
usually of an order of sizep = ab. 0,02, and besides there 
may be additional damping from the device driven by the 
engine. ay will be sufficiently big only if it is due to the 
biggest possible impulse and even then probably only in 2 
special cases, namely a 4-cyl. engine with all cranks in a 
plane (K2 giving subharmonic response at q = 1, whether the 
vibration form is symmetric or not) and a 6-cyl. engine with 
symmetric crank position in triangle (K3 at q = 3/2). The 
direct resonance with this big impulse (although depending on 
another impulse summation, as mentioned) should then be 
situated below half of the maximum revolutions, as the sub- 
harmonic will occur at the double speed, and the system must 
therefore have such an unusually low frequency as can probably 
only happen (in stationary engines) when there is an elastic 
coupling between engine and driven device. In such a case, 
however, the danger of getting these vibrations is completely 
avoided by providing the coupling itself with a certain (not 
very big) inner damping - but of course, it can be important 
to know this, because the undamped subharmonic vibrations 
are of such a violent nature as to make it unbelievable,, that 
they could be completely eliminated by a little damping ‘). 

In ship engines, where the one-node vibration form may 
have a low frequency due to a long intermediate shaft, there 
will in this vibration form be a very big damping from the 
propeller, so that there is normally no danger of subharmonic 
vibrations in these cases. 


A common school example of subharmonic vibrations due 


to impulse distortion is the mathematical pendulum (mass M, 
length L, we = g/L for small vibrations) for which the suspen- 
sion point is given a vertical movement y = f sin2wt, The 
vertical impulse acting on the pendulum mass is then: 


2 
P=- = 4fMw’ sin2wt 
dt 


having a horizontal component Psinx = Px for small angular 
amplitudes x. Disregarding damping (which ought to be negli- 
gible in a pendulum) the differential equation for small vibra- 
tions is 


2 
ML + Mgx = sin2ut- x 
at 


4) In his contribution to litt, 6 the author has given a description 
of such a case. 


— 


and in reduced form 


2 
d(wt) 
with and wi = 
e 


Introducing x = a sin(wt+g) we get a 1. order artificial im- 
pulse of size a4 *€ , free to adjust itself in phase. It is a 


case of subharmonic excitation as the frequency of the vibration 
is half of that of the forcing (here vertical) impulse. For £ = 1 
the artificial impulse will be free to produce increasing vibra- 

tions. For §& somewhat more or less than 1 the artificial im- 

pulse will act as elastic force, giving steady vibrations for 


There will thus be instability (increasing vibrations) in a 
range approximately 
f 


f 
(assuming f << L) 


and outside this range free vibrations will be decreasing, as 
the artificial impulse will here act like a damping. 

The amusing point is, that with the pendulum turned upside 
down the vertical top position may be a stable equilibrium. 
Measuring x out from this position the differential equation 
is now 


2 
2 
d(ut) L 


To get steady vibrations the artificial impulse must here 
necessarily be negative, and we have only one possibility: 


For € smaller than given by this formula the artificial 
impulse is too small to counteract both elastic force and inertia 
force; there will be increasing vibrations and the pendulum will 
quickly fall down. But for € bigger than this value, there is 
stability as the free vibrations will be decreasing. Solving the 
equilibrium condition for the vertical movement f we get 


2 
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and it is seen that the vertical movement required for stability 
is rather big. For &£ = 1 (which means that the frequency of 
the vertical movement is double of the natural frequency of 

the pendulum in ordinary, down-hanging position) we must have 
f > L, and even for the highest possible frequency of the verti- 
cal movement we shall have f > $ L. 


‘ 
| 
i 
j 
q 
4 
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4. Variable mass in one-cylinder and multi-cylinder engines 


In the crank mechanism fig. 4 


the piston movement is with good 


approximation 
m= y = R(i-cos a) +4 R(4-cos2a) 
and = R(sina +5 dsin2 a) 


The inertia moment of the purely 
rotating parts is called m,) and for 


a crank angle a the contribution of 


the oscillating mass 8, t the inertia 
Fig. 4. Crank moment of the crank is: 


mechanism. 


g (sina = g (s in*a+ Asinasinza+4 sin“2a) 


1 2 4,2 4 2 4.2 
(Acosa - cos2a - \cos3a - cos4a) 


We shall disregard the small term an. The constant member 


is included the mean inertia moment m_ = m, 


c 
and the proportion tt. oe is called B The total inertia moment pr. 
c 


cylinder is then 


(14) 
with Bac = B(-A 4, A) for p= (4, 2, 3) 


There is a principal variation twice pr. revolution and two 
smaller variations one and three times pr. rev, Variation compo- 
nents of order higher than 3 may safely be disregarded in practical 
calculations, 


Top 
| | 
\ > 
(4-2 cos pa 
: J. B Pa) 


« 


Besides the mass variation the oscillating mass is pro- 
ducing an outer impulse acting on the crankshaft, and when 
fully considering these two effects active in the vibration 
system, we do not need to think any more of the oscillating 
masses themselves, At steady speed of revolution we, this 
outer impulse pr. cylinder can be calculated (as will be seen 
in the following) as: 


in which latter expression the coefficients 


= 8 RW, (3 for p= (4, 2, 3) 


are the common expressions for the impulses from the oscil- 
lating masses, usually added directly to the cylinder impulse 
Kp of same order, as these big impulses are sinus compo- 
nents with a very small phase angle Yy at the top position. 
Strictly speaking, we should add Kp and Kpo as vectors with 
phase difference Y. We 

Dividing (15) with m, 08 and using q = 7@Q we get the 
impulses from the oscillating mass expressed in terms of 
static amplitudes as 


a (16) 


In multicylinder engines the analysis in free vibration forms, 
explained in section 2, will be correct also when taking the true, 
variable contributions from the oscillating masses into considera- 
tion, These masses are stiffly connected with the constant swing 
masses and do not increase the number of separate masses characte- 
rizing the system. There is still the same number of free vibration 
forms, and they have the same shape as for constant masses, as 
they should (by definition) be calculated for the mean values of the 
variable masses. They have also definite proper frequencies, not 
changing with the size of the vibrations, although (as we shall see) 
slightly higher than for constant masses, 

While in linear multi-mass systems the different free vibra- 
tion components are completely independent of each other, so that 
each of them can exist separately, there may in the variable-mass 


system be a certain, usually very weak coupling between different 
vibration components (this matter is treated in section 10). But in 


} 

: 

q 


each single vibration form A the action of the variable masses back 
upon the vibration component A itself must be the same as the 
action from the variation of the total =m eo this variation being 
imagined equally distributed over all the masses. If it is not so 

in reality, the inequality will have some effect on other vibration 
components, thereby causing the said coupling effect, but not 
changing vibration form A itself. 

If a is the crank angle of cyl. 1, measured from top, its 
variable contribution of order p (p=1, 2 or 3) to the total Im a is 
B cospa, and the of another crank at angle in 
relation to crank 1 is - Bcm™o cos p(a+ 6). It follows that the 
variable contributions from all cylinders are given as p-order vec- 
tors of size B ma” and relative angles p®, giving a resulting va- 
riation coefficient 


p= (47) 


In the equivalent one-mass system representing a single 
vibration form in a multi-cylinder engine the total, variable inertia 
moment will then be 


J=m(i- 2B cos p(a +p)) (14a) 


with m standing for 2m ya and pp meaning the’ phase angle of vector 
8 p when a is still measured from crank 1 in top. 

In the following we shall mainly treat the 2. order component of 
the mass variation and for convenience write B without index in stead 
of B,- Furthermore, we shall omit the deviation angle », now measu- 
ring a out from the "main position" of the shaft, i.e. the position 
at which vector B is at top and the total inertia moment J is mini- 


mum, We have then 
J = m (i - B cos2a ) (14b) 


We can also consider the equivalent system as a one-cylinder 


engine with A= 0 and an oscillating mass 


« 2§< 


otm = Bm (14c) 
Fig. 5 shows the vector summation for the 2. order mass 
variation in vibration form A for the 10-cylinder engine mentioned 
in the introduction, and the whole calculation necessary is shown 
in the following table, in which all inertia moments are given as 
relative values in proportion to the inertia moment of cyl. 1. 


Be was 0,31 for all cylinders, and the resulting B = 0,11. The 


Free vibration forms Crank position 2. order vector 


12345 { 


R=0,69 
th 
56 


Fig. 5. Variable mass in 10-cylinder engine. 


main position, marked I, is at p = 18° from crank 4, which also 


follows from the symmetry of the crank position, 


10 

To propeller 


6 for vibration form A 
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a’ B. 2 ™rel 2 
(2 cyl.) | (pr. cyl.) Ba (pr. cyl.) miei°4 

Cyl. 1+2 1.85 | 1 4.85 

5+6 19 06 2.40 46 

7+8 . 89 .28 4 89 
9+10 | 1,87 58 4 1, 87 
Flywheel | 4.05 0 25 .26 
6.26 

R 0.69 0. 140 


rel 


5. System with evident Corioli's force 


In the first instance one might suppose that the differential 
equation for the now considered variable one-mass system could be 


derived from equation (1) simply by substituting m with m(1- 6 cos2a), 


but this would not be correct because there must necessarily appear 
a Corioli's force acting on the so defined system as an outer force, 


and this must require a new member in the differential equation. 


In whatever way the mass variation is effectuated there must be an 
accompanying Corioli's effect because the inertia of a swing mass 
cannot be altered except by moving the mass to another place in the 
vibrating system. 

This effect will be immediately evident if we look at a con- 


structed, not rotating system fig. 6, consisting of a constant swing 


mes 
2G, Rt 


Fig. 6. Constructed system to show Corioli's force. 


mass m, and an additional mass G, moving in a radial slot and. 
driven (as suggested) by an outer mechanism, independent of the 
vibrations, In order to make comparison with the crank mechanism 
the movement of G,» measured positively outwards from the mean 
radius R, shall be defined as 


cos2wt 


and the total, momentary inertia moment will be 


2 
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2 
J= m, + G(R + y) 


The differential equation for the vibrations in such a system 
is generally (see f.inst. Timoshenko, litt. 8) found by means of 
Lagrange's equation 


d ,OT 
at (Se) 
Using x as parameter we have T = Sax? where J is here 


time-depending only (through y), consequently: 


“Ox Ox 
2 
and Lagrange's equation means that mi¢ in equation (1) should be 
dt 


substituted with 


2 
d(Jx) _ .d“x dx 


in which the second member clearly has the character of the mo- 
mentum from a Corioli's force. It is seen directly from fig. 6 that 
besides the regular inertia forces giving a momentum je there 
must be a Corioli's force 2G sy 


o dt 
as a reaction on the system and therefore positive at the left side 


a giving a momentun? ‘(acting 


of the differential equation): 


dy, dx dJ, dx 
2G (R+y) Gt at’ at 


When f << R we have 


2 = = 
Jez m, + 2G Rf cos2wt = m(4 B cos2w_t) 


with m =m_ + and Bm = 2G Rf 


Expression (18) is then 


2 
d d 
m (4 - B cos2w,t): + 2Bmw sin2w (19) 


and when there is also an outer impulse K sin(wt + Y) acting directly 
on the swing mass, the complete differential equation will be: 


2 
d d d 
m(1-8 cos2w + sin2w t + sx = K sin(ut +y) (20) 


As for the driving mechanism for the moving mass,_ this will 


2 
be acted upon not only by the regular inertia force - cy (which is 
dt 


without interest and could be equalized by a spring), but also by a 
special centrifugal force 
dx \2 
C = (24) 


which may consume work and certainly will do so in the resonance 
case, 

Perhaps the comparison with the crank mechanism would be 
more obvious, if we let G, have the same size as 8, in the crank 
mechanism and let it slide all along a centerline, its distance from 
the center being y = R sinwt, thereby giving a direct impression of 
the oscillating mass g,'8 participation in the inertia moment. The 
differential equation would be exactly the same as (20). But even 
such a system would not, as we shall see, in all respects be equi- 
valent with the crank mechanism, as it is a "freely variable'' system, 
the mass variation being a pure time-function and not influenced by 


the vibrations. Equation (20) is in mathematical terminology a "linear" 


equation with variable coefficients, 


‘ 
4 


- 


6. Correct differential equation for the crank mechanism 


The common crank mechanism fig. 4 is functioning at the same 
time as power drive, pulsator for vibrations and driving mechanism 
for a variable mass, while the two latter functions are separated in 
the system fig. 6. Accordingly, we shall in the correct differential 
equation for the crank mechanism find members not alone correspon- 
ding to the Corioli's force in (19) but also to the reaction from the 
piston rod in its capacity as driving mechanism for a variable mass, 

Using the crank angle a as parameter in Lagrange's equation 
(see also Malkin, litt. 9) we have T = 5J4° with 


J = m(1 - B cos2a) and a= wt +x 
consequently 


and Lagrange's equation means now, that mi in equation (1) should 
dt 
be substituted with 


dt 2 da a ate dt dt 2da ‘dt ate 2 da ‘ dt 
2 
d dx ,2 
= m(1- 8 cos2a) - + Bm sin2a- (22) 


in which the 2. member contains the same Corioli's momentum as in 
(19) (when approximating a with a) although it is difficult to per- 
ceive it as such in the crank mechanism. Besides, this 2. member 
contains a constant times sin2a, which is the regular 2. order 
impulse from the oscillating mass (the same as (15) for p = 2), and 
a term with oS F reminding of the centrifugal force (21). 

The expression (22), meaning the total reaction from the moving 
masses, could also be found directly, without using Lagrange's theory, 


namely by going back to the fundamental equation for the movement of 


the oscillating mass, when A= 0, 


i 
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y = R(4 - cosa) 


and already here introduce a = wit + x. We then get: 


2 2 
dy _ dx dx,2 
at R sina: +i = Rcosa +Rsina ‘at 


The oscillating mass is gy by (14 c) in the general case of an 
equivalent oe: are engine, and the reaction from the inertia 


force - 8s ty on the vibration system is then: 
dt 
a? 4 2 dx 
Rsina =>g'R (sin2 a (w, + (1 - cos2a)- 


which together with the reaction from the constant swing mass, 


a ee and with $8! Rr? = Bm, gives a total momentum on the 


dx 
m (1 - B cos2a)> +Bm sin2a- (w, 
dt 


the same as (22). Having now got a double security in the argumenta- 
tion we can with confidence write the complete differential equation as: 


dx 
m (1 - costa) +B m sin2a- (o, =) K sin(pa +Y) 


(23) 
Remembering that a is not a pure time variable, but = wot + x, 
the equation is indeed very complicated, but we shall see in the 


following, that simple and yet in 1. degree completely correct appro- 


ximations can be made for all practical demands. 


For a full understanding of the action of the various forces 
it will be of some interest also to compare with a rotating 
system as fig. 6, the fixation being substituted with a heavy 
flywheel rotating with constant angular velocity Wo. We shall 
now use a@ = wt+x as parameter, but as J is still time-de- 
pending only, there will be no 0T/a, and we have in expres- 
sion (18) only to substitute dx/dt with da/dt, whereby the dif- 
ferential equation (20) changes into: 


d d d 
m(1-Bcos2w. t): + 2Bmw sin2wt- (0+) + +8x= Ksin(wt + Y) 
(24) 


system: 


30 <- 


thus merely adding a 2. order impulse from the variable and 
now rotating mass. This impulse will be (when transferred to 
the right side of the equation) 


- 28 mw, sinzw t 


while it is only half as big in the crank mechanism with the 
same mass-variation factor 6, namely by (23): 


- sinda 


The reason for the difference is, that in the crank mechanism 
the regular driving of the oscillating mass reduces the "Corioli's 
effect" (if we may call it so in the crank movement), while in 
the freely variable and rotating system fig. 6 the driving mecha- 
nism for G, is independent of the rotation. In order to obtain 
complete similarity between a rotating system fig. 6 and the 
crank mechanism, the moving mass G, should be guided by a 

fixed leading curve as 
suggested in fig. 7, 80 
that the position of G, is 
determined directly by the 
“crank angle’ a=w,t+x, 
ae the movement of G, being 

ALILLALL y = -~fcos2a, The diffe- 


rential equation for such a 
system is completely the 
x. . same as (23), perhaps with 
Lc J exception of the outer im- 
yy pulse K which might still 
be a pure time-function, 
although it is difficult to 
see how it could be trans- 
mitted to a rotating system 


Fig. 7. Rotating variable- ee through the. turning 
angle a. 
——— with fixed It is now seen, that the 


Corioli's momentum: 


gy, da dy da i da 
2G at 2G, Raq (it 26m sinze ar) 


is partly counteracted by the reaction from the roller pressure, 
which is, for f<<R, equal to the total centrifugal force acting 


on G, C= and has a centerdistance c = 2f sin2a, 
so that its reaction (measured positively against the rotation) is 


da da ,2 
- GR(F) 2fsin2a = - Bm sin2a (i 


| 
| 
| 
7 
| 
| 
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which is just the term - at in Lagrange's equation, not present 


when J is independent of the chosen parameter. When J is 
depending on a parameter a, then it must necessarily mean 
that there is a stiff connection between the variable mass and 
the part directly represented by the parameter, consequently 
also a reaction between them, and from this we get the term 
- ah - $5.47, opposed to and half as big as the term Gra 
coming from the first member in Lagrange's equation, apse), and 
representing the Corioli's effect. Thus the different vibration 
systems here described give a nice illustration of Lagrange's 
ingenious theory. 

In the freely variable and rotating system fig. 6 the Corio- 
li's momentum is, by (24): 


dx ‘ 
28 mw, (w, + sin2wt 


while the combined Coricli's momentum and reaction from the 
driving mechanism for the variable mass is, by (23), both for 
the guided system fig. 7 and for the crank mechanism (where 
they cannot be separated): 


2 
Bm(w, + sin2a 


The regular 2. order impulse is in the latter case reduced 
to the half, -B mw,sinZa, but the Corioli's momentum from the 


vibrations is practically the same in both cases, 26 mu, 


times sin2w.t or sin2a, and the vibrations are therefore also 
practically the same, when not specially exited by the 2. order 
impulse. Finally we have in the latter case the member con- 
taining (dx/dt/ analogous with the centrifugal force (21) ac- 
ting independently in the freely variable system, 


dx 

+ 
| 


7. Free vibrations 


To solve equation (23) for free vibrations we shall disregard 
damping and outer impulse, so that the equation is 


2 
m(4 -B cos2a) - + Bm sin2a - (w, + + sx=0 (25) 
dt 


The first question is, whether it will be allowable in the equation to 
substitute a with the pure time variable wt. Seeking a 1. degree 
approximation which will be completely correct for small values of 
the vibration amplitude a, we need only to take regard to all mem- 
bers in (25), which will contain a in 1. power, namely containing 


x, = or = as single factor. Now cosZ2a is in reality, analogous 
dt at 


with the impulse distortion mentioned in section 3: 
cosza = cos2w t - 2xsin2w t 


and as the correction member contains x and should be multiplied 
2 

with ae, we can here disregard it, In the second member in (25) 
dt 


we can first disregard jae ,. next we can disregard the angle cor- 
rection of sin2a in connection with the factor 2u SS, but strictly 
speaking we should not disregard the angle correction in connection 
with the factor w , that means in reality the ordinary impulse 
distortion of the 2, order impulse from the oscillating masses, 
which distortion would be (as an impulse) 


2 
- 2Bmw cos2w.t (26) 


accordingly of 1. power in a. But, as we shall see, this contribution 
will be comparatively small for usual values of q (the number of vi- 
brations pr. revolution), and we shall so far disregard it (but not for- 
get it completely). Finally we can here disregard the main 2. order 
impulse in the same way as we have disregarded all the cylinder im- 
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pulses which will always be present but generally give forced vibra- 
tions only, which have nothing to do with free vibrations, 
The differential equation is now 


2 
m(1-B cos2a): + 28mw sinza- S¥ + 5x = 0 (27) 
at 


in which a may be reckoned as wt. This being so, we shall from 
now on use a = wt as independent variable. Strictly speaking we 
ought perhaps to introduce a new letter for it, but there will be no 
mistakes in using the familar a now as a regular time variable, 
when we remember, that whenever greater accuracy is wanted, we 
must go back to the complete differential equation (23) or (25) and 
here introduce the a-corrections that might be necessary. 

Dividing (27) with mus, and using q = cs we finally get a re- 
duced form 


dx 2 
(1- B cos2a)+ ——5 + 2B sinza +q x = 0 (28) 
da 


analogous with (2), although differing from it by a factor A (q does 
not need to be a whole number in free vibrations), It is remarked, 
that the approximations so far made do not concern the freely variable 
system by (20) or (24) in which system free vibrations are given 
exactly by (28). This may be a guide for the judgement of cases 
where a-corrections will be necessary. On the other hand, equa- 
tion (28) is valid also for a single vibration form in a multi-mass 
system to the same extent as it is valid for a one-mass system, 

Solving equation (28) the next question of approximation concerns 
the P-value which is usually so small that B in 2. power may be 
disregarded in practical calculations. For 6 = 0, that means for a 
linear, constant-mass system representing the mean value of the 
real system, the solution would be an arbitrary sinus vibration 


x = a sin(qa +9) 


2 
Introducing this in (28) the products cos ax and sin2a> 
da 


will deliver new members of order q- 2 and q +2, also being sinus 
terms with phase 9 at the main position, and in the first B-approxi- 


mation the solution must be 


ow 


x = asin(q'a + g)+bsin((q'-2)a + sin((q'+2)a +9) (29) 


in which q' might be different from q, that means, for a given w 
the real number of vibrations pr. rev., q', might be different from 
the given q standing in the differential equation and meaning the 
number of vibrations pr. rev. in a constant-mass system. 
Introducing (29) in (28) we get 3 coefficient equations for the 
‘members of order q', q'-2, and q'+2, all of them still remaining 
sinus terms with phase 9. From the two latter equations b and c 
are found, always containing a factor B®. In their multiplication with 
Bcos2a and BPsin2a they will give q'-order terms (and also terms 
of order q'-4 and q'+4) with factor B°, which may all be disregarded. 
Consequently, we have of order q' only the 2 ordinary members for 
the constant-mass system, so that q' = q. For an arbitrary value 


of a the two other equations then give us: 


b = -k-a and 
30) 
-2 4 q(q+2) ( 
with k = and 1 = 


and the vibration consists of a major harmonic of order q and arbi- 
trary amplitude a, and smaller harmonics of order q-2 and qt2 
and amplitudes (30) proportional to a, All of them have the same 
phase @ at the main position, a = 0, but 9 is itself arbitrary (and 
also shifting for each passagé of the main position, if q is not a 
whole number). As the vectors representing the 3 components are 
rotating at different speeds, their relative angles are continually 
shifting, but still they are firmly bound together by their com- 
mon position when the crankshaft is at the main position, which 
determines the timing of the variable-mass effect. 

For the 3. order component of the mass variation (with coef- 
ficient B, by (17)) a similar calculation is easily made. The dif- 
ferential equation will be similar to (28) only with 3a in stead of 
2a, and also coefficient 3 in stead of 2 in the second member. 
There will be secondary vibration components of order q-3 and q + 3 


and size: 


| 
| 

| 

| | 
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b=-a- +p, c = (30a) 


“2 
In the present case with a 2, order mass variation B = 0,114 
and q = 8 we get 


10 


k= 0,094 1=0,44° = 0,122 


Fig. 8 shows the 3 harmonics and the resulting vibration during 
4 revolution in the case p= 0, all 3 components being sinus curves 
with O-point at the main position I. For the sake of clearness, how- 
ever, the secondary components are much exaggerated and correspond 
to B = 0.4. The resultant (which would be very similar for other 
values of @) consists alternately of two shorter and two longer vibra- 
tions, which was also clearly seen on the torsiograms, corresponding 
to a ''momentary frequency" (understood rather as momentary phase 


velocity) : 


m "(4-8 cos 


n (4 +48 cos 2a) (34) 


varying in the real case 5.5% above and below the mean value. 


I I 
| 


Fig. 8. Free vibration with variable frequency (B = 0.40). 
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It might be imagined, that eq. (341) could be used as basis 
for a very simple and yet correct vibration calculation, namely 
d(wt) 

dt 
in the expression x = asinwt) is varying as (31). In our notations 
we would write 


by expressing that the phase velocity of free vibrations ( 


x = asing® 


where ® is a parameter given by 


= i+ cosza 


The free vibration would then be (for p= 0) 
x = asing(a +48 sin2a)= a(singa +48 qsinza- cos qa) 


= a(singa - $84 sin(q - 2)a +28q sin(q + 2)a) (32) 


So this simple calculation, altogether filling 3 lines, comes 
rather close to the solution (29) and (30) found above through 
lengthy operations, and in some practical applications (32) may 
be sufficient. But the method is correct only under the pressuppo- 
sition of a high q-value, that means many vibrations pr. mass 
variation period, and for usual q-values the difference between 
the coefficients by (30) and (32) is not unimportant. Above all, 
we should not from this starting point be able to proceed further 
with the theoretical investigation as done in the following. 


Strangely enough, (32) would be correct, if the Corioli's 
force were half as big as it is, which is easily seen from (30) 
where the ciffer 2 in the numerators stand for Corioli's force. 


| 
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8. Fictive forces. 


Subharmonic vibrations due to mass-variation 


The procedure followed above in solving equation (28) is, 
mathematically seen, a simple iteration depending on B <1 (as it 
must necessarily be in our case), The main vibration is the first 
approximation, the secondary vibrations the 2. approximation, and 
in this way we could continue in rapid convergence as f is usually <<1. 

This procedure can be interpreted physically in a way very use- 
ful for the further investigations. When introducing the main vibration 
x = asin(qa +) in the differential equation the new members formed 
of order q-2 and q+2 (that means the members not belonging to the 
linear parts of the differential equation) can be considered as "fictive 
forces" in the distinct meaning, that when we add these forces (and 
their possible later derivations in the following approximations) on 
the system as outer impulses, the system is to be calculated as a 
linear system with constant swing masses equal to their real mean 
values. It is in the same way that the name "fictive forces" is 
used in relative movement. They are in fact very real forces acting 
on the given, real system, but ''fictive'’ when we apply them on a 
simplified modification of the system in such a way as to make the 
result correct in the real system, 

In our case the main vibration or any forced vibration (calcu- 
lated in the system as linear) is a first approximation for the 
solution, and its fictive forces represent, in their action on the 
system as linear, the next and usually sufficient approximation. : 
Looking at the matter in this way we get a clear understanding of 
the various vibration phenomena, and also a great simplification 
of all detail calculations, partly because the simple expression for 
a fictive force is a kind of "operator" to be used many times in 


actual calculations, partly because these operators are vectors of 


given orders and phase angles which are easy to handle, and we 
get all the advantages of vector calculation. Especially for multi- 


mass systems (and all vibration systems of real practical importance 
are multi-mass systems) it would seem impossible to get through the 
calculations in any way than this. 

Dividing equation (28) with * its 3 members - an inertia force, 
a Corioli's force and an elastic force - will be expressed in terms 
of "static amplitudes''. We wish to consider the new members, the 
fictive forces, as outer impulses and must therefore transfer them 
to the right side, changing their sign. From the q-order main vibra- 
tion we get then, by the few members of order q-2, a fictive force 


of this order, expressed as static amplitude: 


5 
and a fictive force of order q+2: - (33) 
1, 1, 


The fictive forces shall here always be denoted by a dot over 
the letter denoting the vibration which has produced them, and with 
index denoting the order of the fictive force itself. (33) is valid 
not alone for a one-mass system but also for a single vibration 
form in a multi-mass system with the appropriate value of B, be- 
cause the action of the fictive forces from the different cylinders 
back upon the same vibration form is depending on the same kind of 
vector summation (13) as is also used in calculating fB by (17). 
(Further about this and also about their action on other vibration 
forms, see section 10). 

The fictive forces (33) are both sinus components with the 
same phase angle 9 at the main position as the main vibration a 
which has produced them. They will themselves produce ordinary 
forced vibrations in the system, now considered as linear. But 
free vibrations mean, that the q-order main vibration is at reso- 
nance, consequently a4-2 is at a relative frequency € = 5-2 


q 
giving a magnification factor 


2 

4 
Fs = 
1- 4(q-1 


la . 

: 

| 

| 
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and a forced vibration 


-2 
q- 


in accordance with (30). Similary we get for the q+2 order vibration: 
cs S28. (- ) = a- 


Thus the secondary components of the free vibration are in 
reality forced vibrations coming from the fictive forces which are 
produced by the main vibration. Therefore the secondary compo- 
nents have also the same phase as the fictive forces (and as the 
main vibration in the main position when q > 2), and they have dif- 
ferent sign because one is forced in under-resonance and the other 
in over-resonance. 

Fictive forces are not produced exclagively by resonance vibra- 
tions, but by any forced vibration which may exist in the system. 

F. inst., if we have a forced vibration of order q-2, x= bsin(q-2)a, 
caused by an impulse of this order, we find by introduction in (28) 


a fictive force of order q 


2 
= - + 2(g-2) a2 
(34) 
and similarly from a c-vibration of order q+2: 


accordingly the same "coupling factors" as by (33) from a q-order 
(free or forced) vibration to its accompanying fictive forces of order 
q-2 and q+2. It is very satisfactory that we have thus the same 
coupling factor going both ways, but this would not have been the 
case if we had not got the member expressing Corioli's force in 
the differential equation. 

If b and c in (34) are themselves just the secondary compo- 
nents caused by the main component asin(qa +)) of a free vibration 


through its fictive forces (33), the secondary fictive forces (34) will 


q-2 8 

| 4 
| : 
| : 
| 3 
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be a kind of "reflection" back upon the free vibration. They will be 
not alone of the same order but also the same phase angle as this, 
and they will therefore act like elastic forces, changing the resonance 


frequency a little, Their size is together 


+2 1 +2 


q 


q- 


accordingly an outer impulse opposed to the vibration amplitude, there- 
fore adding to the elastic force and increasing the frequency. Or we 
can say, that the order number of free vibrations is indeed increased 
to q' as anticipated in (29), but the difference depends on 6 in 2. 
power and was therefore disregarded in first approximation. If we 
now introduce (35) as outer impulse in the differential equation (28) 


linearized, we get 


2 
= -a- 3p" = (4 +387) 
q 


and the mean frequency of free vibrations will be 
~ 
(4 +7} B ) (36) 


The same result is found from (31) if we add the next member 
in the series 


= (145 8*cos“2a)= (145 B*(4+cos4a)) 


giving a mean value (36). Also this calculation presupposes a rather 
high q-value, otherwise the effect is not fully developed, and as seen 
from (35) the correction is in reality sligthly less than given by the 
last figure. 

In the present case with B = 0.11 the frequency increase is only 
0.2 % and without practical importance, but it is of some interest to 
note that also in the variable-mass system the mean frequencies of 


free vibrations in the various vibration forms are absolutely constant, 


| 


« 


not changing with the size of the vibrations as in other kinds of 
(true) non-linear systems. (But the introduction of fictive forces is 
just as useful in these cases, see end of section 9). 

For the other components of the mass variation we shall, of 
course, get fictive forces similar to those calculated above and 
acting as additional outer impulses. It follows, that the different 
components are in no way interfering with each other and can be 


calculated separately as the case may be. 


By means of fictive forces we can easily solve the problem 
of subharmonic vibrations due to the variable-mass effect, 
which is the sole subject of the before mentioned highly 
theoretical papers (litt. 1-4), although the authors do not 
clearly mention that it is in fact subharmonic vibrations they 
calculate. They only speak about a "'splitting-up" of the re- 
sonance points in 2 neighbouring frequenciés at which steady, 
free vibrations can exist, and between them a range of "'in- 
stability", as already described in sect. 3 in connection with 
subharmonic resonance vibrations exited by a cylinder im- 
pulse of the double frequency, but they do not clearly state, 
that in the present case (where cylinder impulses are not 
necessary) the vibrations are directly exited by the mass va- 
riation of the double frequency through a fictive force produ- 
ced by the vibration itself. Therefore the energy supply to 
the vibrations remains a mystery to them, because, how can 
instable vibrations arise and increase without impulses sup- 
plying the energy ? Of course the mass variation is always 
in rotating systems accompanied by an outer impulse, but it 
is here of double frequency of the vibrations and cannot sup- 
ply energy to them (except through impulse distortion which 
these authors do not consider, and besides we can have sub- 
harmonic excitation even if the total outer impulse is 0). 

We shall here first treat the matter for a freely vari- 
able system for which the free vibrations are exactly given 
by eq. (28) with a as a pure time variable (whether the 
system is rotating or not does not matter, because the outer 
impulse from the mass variation will here also be a pure . 
time function, not subject to impulse distortion). We can 
without trouble generalize eq. (28) to be valid for any of 
the 3 mass-variation components of order p = i-3, namely 
in the form: 


2 
(1-B ,cospa) - + + = 0 (37) 


The subharmonic condition is p=2q. Introducing 
x = a sin(qa + 9) we get a fictive force, which expressed in 
full will be: 


x 
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=a: (38, - sin(ga - @) = sin(qa - 9) (10a) 


This force is in every respect analogous with (10). It is 
in resonance with the vibration, and the vibration is (due to 
the conversion of the sign of p) free to adjust itself to any 
suitable phase difference. If there is damping, we have the 
limit condition 


p= 58 (14 a) 


and if p < 38, the range of subharmonic vibrations is given as 


- (12 a) 


For p = 0 we have €£ = 1+ 56, which is the main result 


of the said papers (litt. 1-4). In practice this result is quite 
insufficient, partly because we have always a damping, and 
partly because we have always a big cylinder impulse of the 
same order as the mass variation and subject to impulse 
distortion, producing the artificial impulse (10). Of course 
there is great difference between "artificial impulses" due to 
impulse distortion and "fictive forces'' due to mass variation, 
and we shall here retain these different names, but there is 
also a great similarity between them, particularly the feature 
that their size and phase in multi-cylinder engines are de- 
pending on the same kind of vector summation, namely ay by 
(13) and B, by (17), In the formulae we can therefore without 
mistake use the same kind of notation for them, as they are 
in any case to be added directly if they are of the same 
order. 

Thus, in calculating subharmonic vibrations in crankshafts 
we shall add the artificial impulse (10), -qapy. a, and the 
fictive force (10a), “3 Bp: a, if they are both of a size to be 
considered. a, may exist for any order number p, while Bp 
is here considered only for p = 1-3. But in the cases where 
the mass variation is considerable we have at the same time 
a resulting oscillating-mass impulse, also subject to impulse 
distortion. It should therefore be included in apy, what it 
also is if we in calculating apy by (13) already include the 
oscillating-mass impulse pr. cyl. in the cylinder impulse, 
As there is a certain resemblance between the two effects of 
the variable masses (one through the inertia- and Corioli's 
force, and the other through impulse distortion) we need 
clear distinctions, and we shall here provide thatapy always 
means the resultant of cylinder impulses only, while we shall 
use the notation a%, for the resultant of cylinder impulses and 
oscillating-mass impulses, These latter are given by (16) for 
a one-cylinder engine, and the expression is valid also for 
multi-cylinder engines with Bp in stead of Ppc. Furthermore, 
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the phase difference between cylinder impulse and oscillating- 
mass impulse is negligible, 


so that we have generally: 


a* =a - 
pu pu (16 a) 
In the subharmonic case p = 2q we have 
a* sa - 4 6 
pu pu q 
and the total fictive force will be 
a' = - a+ (qa* ) = - a-+(qa ) (10b) 
q pu pu 


and for the 


on 
ne 


our results 


that means, 
the cylinder impulses Kp after adding only half of the oscilla- 
ting-mass impulse, a curious but quite satisfactory result. 
We have now for the limit damping: 


specially for p = 2 


the force is the same as would be calculated from 


1 
1 
om (14 b) 
1 
range of the vibrations: 
4 


This formula includes all possibilities p = 2q and combines 


from sect. 3 with the results from the previous 


papers on the subject (litt. 1-4), which do not consider im- 
pulse distortion and damping, but agree on the member 4 


In litt. 5 (from a much later date) due regard is taken to 
damping but not to impulse distortion. 


Regarding the question of energy supply it is now obvious, 
that in the freely variable system (whether rotating or not ) 
the energy must be supplied by the driving mechanism, for the 
moving mass, 
at resonance it performs a work pr. q vibrations 


In these cases the fictive force is -a- Bp, and 


4 2 2 
A= 54956, a 


J 

ee 


+ « 


The driving mechanism is pr. p periods performing a 
work against the centrifugal force (21): 


where y = -fcospa ae = f pw sinpa 


and 2G Rf = giving 


For the crankshaft system the matter is more complicated. 
We have here the total fictive force (10b). The first part of it, 
-a-+qafu, is the primary artificial impulse which, as shown in 
sect, 3, is "braked" by its accompanying secondary artificial 
impulse. The other part of it, - a: 4Bp, is the same as in the 
freely variable system and gives the same impulse work A. 
Braking work to the same amount must then be supplied by the 
members in the correct differential equation for free vibrations 
(25), which are not effective in the freely variable system. We 
shall, however, leave out the variable-mass impulse which has 
been included in afy. In generalized form eq. (25) is then: 


2 
(1-B ,cospa) + PP ,sinpa +3 PB sinpa + =0 (38) 


where a in cospa and sinpa in reality means %ot+x. Retaining 
the letter a as a time variable, cospa and sinpa should be, 
respectively 


cos pa - pxsin pa and sin pa + px cos pa 
We have then 3 new members: 
N = pf sinpa +p 28 cospa sinpa (39) 


one from the inertia force, one from the Corioli's force, and 
one from the "centrifugal force", which all contain the product 
of two of the functions x and its derivatives ang therefore 
will give a braking momentum proportional to 


p~ 
; 2 
1 2.2 4 22 
A= -5pG [singe + (4-sinpa)da = 
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For x = asin(qa +4) we have, similarly as above: 


2 
da 
dx ,2 422 
(=>) =5aq (1 - sinpa) 


In order to get dimension of static amplitudes eq. (38) 
and N should be divided with q*, and we get the braking work 
pr. revolution: 


2n 


The "centrifugal force" gives here a braking work -A, that 
means a forward pull to this amount, but the inertia- and 
Corioli's forces give together a braking work 2A. It follows, 
of course, from the energy principle, that in crankshaft vibra- 
tions any kind of impulse work, whether from real or fictive 
impulses, must be counterbalanced by a braking work, and 
there is nothing surprising in this, but we have found it 
advisable in this case to show the calculation in details as a 
final control on the above developed formulae, especially 
formula (30) (and the corresponding (41) in next section) which 
are fundamental for practical calculations, Whether such ex- 
pressions are a little more or less complicated (f.inst. whether 
the coefficients should be taken by (30) or by (32)) does not 
mean anything in the calculation work, but it is very satisfac- 
tory - when we have to be content with a 1. degree approxima- 
tion - that this approximation is at least completely correct 
within its own limits, 


We can now also make up for the member (26) which we 
inadvertently disregarded in the calculation of free vibrations, 
sect. 7, and which means the impulse distortion of the os- 
cillating-mass impulse. But in real engines, even in "free" 
vibrations, we cannot have such a member without also having 
distortion of the cylinder impulse of the same order. The 2 
impulses together are given in (16a), and they produce a 
fictive force a-4 pat, which should be added to the main fic- 
tive forces (33). Expressing these in a more general form we 
finally get: 


‘9 


The expression is valid for all values of q and p. Only 
for p = 1-3 we have a mass-~-variation effect, and for p > 3 
the first member in the parenthesis disappears, but the other 
member remains as pa,,, which is the usual impulse 
distortion of a p-order impulse by a q-order vibration. For 
p = 2q the lower fictive force (33a) reduces to the subharmo- 
nic value (10b), although with opposite sign because (33a) 
gives al, = -aj. In the subharmonic case we need the correct 
expressidn (10b) because the two forces will probably be of 
the same order of size, a8 apy will be big due to the extra- 
ordinarily low frequency in such cases. But in normal engines 
running at much higher q-values a},, will usually be small in 
comparison with B,, and in the general theory we shall there- 
se use expression (33) or generalized as the first member in 
33a). 

(33a) shows in any case, that variable-mass effect and 
impulse distortion are two things which can be considered 
separately - except in the subharmonic case, where their com- 
bined effect is decisive for the existence or non-existence of 
the vibrations. 


9. Secondary resonance, 


The use of fictive forces in true non-linear systems 


It is obvious, that "non-linear" vibrations of the kind (29) and 
fig. 8, having in their free state secondary components of order 
q-2 and q+2, can be excited by impulses of one of these orders, 
when the main vibration of order q is at or near resonance. All 
three components are present in the free vibration, and when im- 
pulse K, is in resonance with the main component, impulses K,-2 
and Ko+2 can be said to be in "secondary resonance" with the 
vibration, as they are obviously able to deliver impulse work on 
the total vibration although in much smaller degree than the im- 
pulse K, in "direct resonance", 

The case most likely to be met with in practice is when the 
impulses K, in a multi-cylinder engine give a very small resulting 


static amplitude a so that this resonance point would be consi- 


dered unimportant ad the usual calculations, while the impulses Ko2 
give a very big @(g-2)s° the direct resonance of which would be 
situated well above maximum revolutions and therefore also disre- 
garded in usual calculations. If in such a case the system has 
also a considerable 2. order resulting mass variation 6, the vi- 
bration can be excited to a rather big size by the impulse Ko-2 
in secondary resonance. This was just the case mentioned in the 
introduction and calculated for mass variation in section 4 and 
for free vibrations in section 7. 
The excitation in secondary resonance is very easy to cal- 
culate by means of the theory of fictive forces. The impulse Ko -2 
with static amplitude (q-2)s will at the resonance point w, =~ 


give a forced vibration of order q-2 and amplitude: 


2 


which according to (34) gives a fictive force of order q: 


« 


i.e. the same transformation factor, -k, as between the components 
of the free vibration given by (30). The main vibration is then given 
as 


(1- + 


with suitable phase difference between a and b!, as the fictive force 
is acting completely like an ordinary impulse, and the vibration is 
free to adjust itself to the right phase difference, Therefore the 
whole stress curve about the resonance point is of the same shape 
as in linear vibrations, 

The matter seems very simple, especially as the transforma- 
tion factor (41) is equal to the relation between the components of 
the free vibration, and yet it would not be quite correct to think, 
that the impulse Ko-2 is acting directly on the secondary compo- 
nent of the vibration which is, after all, a rather loose and just 
"secondary'' apparation, The real thing is, that the impulse K,-2 
is in the first place producing a very definite and constant forced 
vibration of its own order, and it can then be imagined (as the 
calculations also show), that this forced vibration through the inner 
mechanism of the system can produce exciting forces of other or- 
ders. The secondary component of the free vibration is only an 
outward sign, that such a mechanism exists, which will also work 
the other way. 

Therefore the forced vibration by has nothing to do with the 
b-component of the excited vibration. by is fixed at the same phase 
as the resulting impulse 4(q-2)s’ while the b-component changes 
phase together with the main vibration as the degree of resonance 
may require. It is then so, that the exciting force bg is depen- 
ding solely on by which is constant in the interval to be consi- 
dered (or only slightly increasing with the revolutions) and does 
not require any impulse work. But at resonance, when the main 
vibration a follows in phase $ after b' ,the b-component will fol- 


low $ after (g-2)8’ and the impulse can then deliver work on it 


it 


- 49 - 


(although not the full impulse work consumed by the vibrations, as 
shown below). 

Usually the impulse K, in direct resonance cannot be com- 
pletely disregarded even if its resulting static amplitude a__ is 


qs 
comparatively small. 60 should be added geometrically to bo 


and this is easily done by considering the crankshaft in main po- 
sition. Here the fictive force is always in phase opposite to the 
forced vibration which has produced it. In practical calculations it 
will be more convenient to use the ordinary vector diagrams for 
impulse summation drawn for crank 1 in top. Then, if » is the 
angle by which the crankshaft in that position has passed the main 
position, bg will be at an angle 2, past the opposite of a 


q-2)s 
(see also fig. 10). 


When using the corrected value for the coupling factor 
(33a), which is also applicable to (34) because the coupling 
is the same both ways, and further introducing a general 
expression for the forced vibration b4we get a complete 
formula for fictive forces of order q, namely 


2 + 
1, 9-P 14 
b! (29 * 44 


From this it is seen, that a fictive exciting force for a 
q-order resonance vibration may be produced by any 2 im- 
pulses for which the sum or the difference of the order num- 
bers is equal to q. The conditions are, that one of the im- 
pulses should give a big 4g-summation and consequently a 
comparatively big forced vibration, while the other impulse 
should give a big ay -summation, consequently a considerable 
impulse distortion or mass variation to be worked upon by the 
forced vibration. 

Most of these cases are, however, subtleties, and in 
practice we shall probably never meet this kind of excitation 
of any important size except in the case p=2 and impulse 
of order q-2. Even in that case the conditions are narrow, 
namely that both B and a q-2)s should be big, and the effect 
depends on the product of is: 

But then we have also the subharmonic cases p= 2q 
which, seen from this point of view, are curious singula- 
rities at which the resonance vibration itself plays the role 
of the otherwise necessary forced vibration, although with 
negative sign, namely for 


= 
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When the vibration comes up at all, its disturbing effect 
is much bigger than that of an ordinary forced vibration, and 
there is now only one condition to be fulfilled, namely that 
apu Should be big (correctly qapy -4 Bp). This requires a 
rather low p-value, but not necsessarily as low as 2, We 
shall therefore not give any definite rule for these cases but 
only say, that any case p=2q at a comparatively low p-value 
should be regarded with suspicion, untill it has been calculated 
at least approximately. 


We shall finally investigate the energy balance in the case 
of excitation by an impulse Ky-2. Assuming y = 0 (its value 
is without influence on the size of the vibrations) and writing 
ag Short for a(q-2)s we have at sharp resonance: 


Impulse sin(q-2)a 
2 
Forced vibration by’ sin(q-2)a with b, = a, Weiy 
Fictive f -ka.- sina with = 
ictive force a, a 
Main vibration cosqa 
b-component -k a+ cos(q-2)a 


The work done by the fictive force on the main vibration is 
pr. revolution (in relative measure): 


qxkaa, 


while the work done by the real impulse on the b-component 
is only 


(q-2) xkaa, 


The missing work 
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is in the freely variable system delivered by the driving mecha- 
nism for the moving mass. In (dx/dt/ determining the centri- 
fugal power (21) the double product of the velocities of the 
main vibration and the forced vibration will namely give a 
member with sinZa which by multiplication with the sin2a in 
dy/dt gives a constant of the right size. 

In the crankshaft system we have the 3 "braking members" 
(39), for p = 2: 


ds dx dx \2 
N= 28x sin2Za + 4Bx—— cos2a + sin2a 
dae da da 
Here the products of the main vibration and the forced 
vibration give us the following members of 2. order: 


from : (q- 2- q) cos2a =- ab, cos2a 


from : ab, q(q-2) sin2a 
2n 


The braking work {3 da is now from the inertia- and 
q 
° 


Corioli's force together: 


ab 
4 (2 _ 
ZnB (q_-2q+2-2) = 2nafa. 


and from the "centrifugal force": 


» 
a 
" 
N 
> 


ab 
2x8 a(a-2) = ~2nkaa, = -A 
q 


so that we have the same as in the subharmonic case, sect. 8: 
The centrifugal force gives a forward pull, or braking work 
-A, while the two other forces together give a braking work 
2A. The total braking work A from the 3 members is pro- 
portional to a, because it is proportional to bi, and it is 
transferred to impulse work without the direct action of ag . 
But ag is delivering the rest of the impulse work, thereby 
consuming a braking work to the same amount through usual 
impulse distortion. 


| 


In true non-linear systems in which the variable coefficients 


of the differential equation are depending solely on the vibration am- 
plitude the use of fictive forces will give the same simplification 

of the calculations as in the freely variable (in mathematical sense 
"linear'') and the guided variable systems, We shall here just 

show the application of fictive forces on the simplest possible non- 


linear system, namely a one-mass system with variable stiffness 


s = 8 (1 + 6+ x) where |o-x|«4 
2 8 
Using a = wt as independent variable with w, = = the dif- 
ferential equation for free vibrations is; 


2 
— > + x(1+6x) = 0 or x'+x=-0+x 
da 


so that the fictive force is simply - éx*. 
For the free main vibration x= asina we have fictive forces 
of order 0 and 2, namely 


s -4 ba" and a = 56a’. cos2a 


consequently there is a 2. order component in the free vibration 
besides a constant displacement of less practical importance. Se- 
condary resonance and subharmonic excitation is here the same 
thing, and it would be impossible to say which name is the more 
correct. In the freely variable system secondary resonance means, 
that the exciting force is produced by a forced vibration, while 
subharmonic excitation means, that it is produced by the main vi- 
bration itself, but in the true non-linear system the exciting force 
is produced by the product of the main vibration and a forced vi- 
bration. This latter, forced by a 2. order impulse ar5° is 


and the total vibration is (disregarding the constant displacement): 


x = asin(a+) + sinza 


4 4 
b, = ‘a = a 
1-28 28 3 “2s 


« §3 
giving a 1. order fictive force 


with free phase difference in relation to a, so that the equilibrium 


condition in the resonance case is 


p =Z6a,, (14 c) 


It will not be possible to obtain this result in an easier way. 

So in this method we cannot distinguish between freely variable and 
true non-linear systems, because in both cases the fictive forces 
represent everything besides the "linear" forces coming from the 
constant members in the differential equation. But the fictive forces 
can be of very different character in the different cases, 

The fictive force (10c) has the character of subharmonic excita- 
tion as it is proportional to a like the force (10a) in a freely va- 
riable system, but it has also the character of secondary resonance, 
as it is proportional to an a. like the force (41). Therefore, the 
presence of a forced vibration and the absence of damping above a 
certain value (11c) in relation to this are the conditions for the vi- 
brations, and although the exciting force itself is only appearing at 
the same time as the vibration, this will inevitably arise at the 
slightest, unavoidable disturbance, 

In the case s = (4 +e° x”), for which the differential equation 
is called Duffing's equation, we get similarly: 


b4= e+ and p= (14d) 


for subharmonic vibrations caused by a 3. order impulse a36° 

In this case a given (small) damping can prevent small vibra- 
tions but not bigger. The vibrations will not come up by themselves, 
but they can be started by a single push, and they will then increase 
rapidly. 

All of the vibration energy is in these cases supplied directly 
by the impulse, as there is no special mechanism involved in the 


variation, 


4 
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In relaxation vibrations following van der Pol's equation 
2 
x' = e€(1-x)-x' +x = 0 


we have a total fictive force e(4-x°) x', Introducing 
= asin(a+p) we get a 1. order fictive force 


a cos(a (1 - +o) ) 


2 2 
= acos(at+@)- (4 +e a> cos(a + cos2(a +m) 


2 2 2 
= acos(a+@)- (1- =) +e a cos(a +9)> ze cos(a + 


The vibration is then creating its own impulse through some 
special mechanism, also supplying the energy from outside, just 
as in the case of subharmonics in the crankshaft system. This 
artificial impulse is of constant phase difference, for a < 2 
always leading the vibration by 1/4 period and causing increasing 
vibrations until a = 2, at which point the 1. order impulse is 0, 
and steady, free vibrations can exist, although distorted by se- 
condary components of forced vibrations caused by the fictive 
forces of higher orders. For big values of ¢« these distortions 
will be big and also reduce the frequency of free vibrations, 
but the maximum amplitude remains at a = 2, (see litt. 11, 

p. 4152 f). 


Also for calculation of true non-linear multi-mass systems the 
method of fictive forces will be useful. We shall not treat the matter 
here but only say, that in stead of the summation 2m a repre- 
senting kinetic energy comes here the summation 2 s° a? represen- 
ting potential energy, s meaning the stiffness of the various parts 
of the system and d their difference amplitudes, By this it is easy 
to determine an equivalent one-mass system for each vibration form, 
although the matter is quite different from (or nearly opposite to) 
the problems of freely variable systems, as the fictive forces in the 
true non-linear system are not single, free forces but difference 
forces between neighbouring masses. Therefore we must here in 


stead of the usual rope polygon fig. 3, representing amplitudes over 


the elasticities as abscissae, use the complementary rope polygon 


§§ « 


representing forces over the masses as abscissae. The fictive 
"forces" (force differences) will all be in same phase, so that we 
have no trouble with phase angles in this case. On the other hand, 
they will not be proportional to the "vibration form" (force form) 
which has produced them, but rather to some higher power of it, 
and they -will therefore contain components also in the other vibra- 


tion forms causing coupling effects similar to those treated in the 


next section. 


40. Coupled vibrations in different vibration forms. 


Coupled subharmonic vibrations. 


The vector diagram fig. 5 can now be considered as repre- 
senting fictive forces of order q-2 or qt+2 in their action back 
upon the vibration form of the q-order vibration which has pro- 
duced them. The full understanding of this requires, however, a 
little reflection. 

The fictive forces from each individual cylinder has not, like 
ordinary cylinder impulses, a given phase when the crank for that 
cylinder is in top position, but they have in this position the same 
phase as the vibration which produces them (correctly: they are in 
phase opposite to the vibration as they have opposite sign). From 
the time when a crank is in top until the time for which the vector 
diagram is drawn the crank has turned an angle @, the vibration 
vector has turned an angle q®, and the fictive force vector an 
angle (q + 2)@, the difference being + 2. At that time the vi- 
bration vector is in some phase which in any case is common for 
all cylinders, and the diagram gives therefore the correct relative 
position of the fictive force vectors, Their size is proportional 
to Bo 4.) and as the vectors are set up in the diagram as p.A., 
they represent impulse work vectors, and the resultant R repre- 
sents the resulting fictive force expressed in terms of static am- 
plitude as calculated in the formulae above. 

There is now only one thing more we need to know, namely 
the real phase angle of the resulting fictive force at a given time, 
but as fictive forces and mass variations are composed by the 
same vector diagram, the above rule for a single cylinder can be 
extended to the whole vibration form, so that the rule is: When 
the crankshaft is in the main position, the resulting fictive force 


vector is always opposite to the vector of the vibration which has 


produced it, Although the two vectors are rotating at different 


speeds, they are timed together at the main position, at which re- 
sultant R, fig. 5, has phase 0, This means, that if the vibration 


«it « 


has here the phase 0, the resulting fictive force has phase x, and 
if the vibration has another phase, the fictive force is always oppo- 
site to it. The argument is valid both for fictive forces of order 
q-2 and qt2. 

So far there is nothing new in this, but it is now seen, that 
we can easily calculate the action of this set of fictive forces on 
another vibration form B, namely by using the same vector dia- 
gram fig. 5, only with vectors of size BAAS, in stead of Ban ° 
We shall then get a coefficient of "coupled mass variation", ana- 
logous with (17) but here considered only for the 2. order mass va- 


riation: 


= (6 c™4,4,) 


(43) 
ab may 


from which we can calculate the fictive force * in vibration 


form B produced by vibration a, in vibration 2 ban analogous 
with (33): 
2 
1 
The notation is a little troublesome, but the idea is simple 
enough. Bab by (43) means the action from vane form A to vi- 
bration form B, and we have therefore put zm) in the denominator, 
but in order to get the fictive force (44) expressed in terms of static 
amplitude in that vibration form we must further multiply with the 
square of the relation between the frequencies of the two vibration 
forms. 
In practice it is more the opposite way which is of interest,. 
as we wish to know all fictive forces of order q in the vibration 
form A, whether coming from this or from other vibration forms. 
We need here the mass variation coefficient "from B to A": 
2 (B.m_4,4,) 


ba 2 


(45) 


which will give an additional fictive force, similar to (41): 


| 


%- 


+2 
(46) 
where d,means the real, forced vibration of order q + 2 calculated 
separately in vibration form B and therefore not to be expressed 
by a general formula like (40). 9(q4 2)s will be smaller in the 
higher vibration form and the magnification factor closer to unity. 
The phase of qo is found by the usual rule, opposite to the forced 
vibration at the main position which, however, will be different 
from the main position for ® (or Dad as it should now be called). 
Such coupled vibrations will probably never be of importance 
in crankshaft vibrations except in very special cases, but it is one 
of our purposes here to anticipate all special possibilities. As for 
coupled vibrations these possibilities are easily discovered by 
drawing the simple vector diagrams for the P-values (45). They 
could also be found by calculating the exact vibration form for va- 
rious positions of the shaft and then note the variation in the shape of 
the vibration form, but this would, of course, be very troublesome. 
Fig. 9 shows the working 
of the coupling in a case where 
| there is a high Boat The full 
drawn curves are the real vibra- 
tion forms of free vibrations, 
when the crankshaft is (stationary) 


at the main position for Puce and 


at a position 90° from this. 
During rotation there will be a 

Fig. 9. The working of continuus shift between the two 

coupled vibrations. modifications of the vibration form, 
and the free vibration will contain 

a q-2 order component of vibration form B, which can be considered 
(and calculated) as a forced vibration produced by the force (44). 
Oppositely, a q-2 order impulse giving a big a, in vibration form B 
would be in secondary resonance with the q-order vibration in A, 
Boa may even be 0, and we have then no frequency variation, only 
alternately increasing and decreasing vibrations in the two cylinder 


groups. 


| 
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We can also get coupled subharrnonic vibrations in two 
vibration forms A and B, caused by the 2, order mass vari- 
ation (or the cylinder impulse Kz), namely at a speed of 
revolution: 


wo, to, 


when wa and wh correspond to the proper frequencies of the 
two vibration forms and these between them have a conside- 
rable Bab The numbers of free vibrations pr. revolution in 
the two vibration forms will namely be: 


Zw, 
a wo, +o, b w to, 
q, (a, 2) 
so that + = 2 or 
(qa, - 2) 


which means that each of the two vibration forms will produce 
fictive forces in resonance with the other vibration form. 
Furthermore, due to the conversion of the sign of the q-va- 
lues there is also a change of sign of the phase angle 9 when 
going over from the vibration to its accompanying fictive forces, 
and consequently there is free possibility of phase adjustment 
in the neighbourhood of resonance, 

When the damping is small enough we have then increasing 
vibrations in both vibration forms at the same time as each 
component will excite the other in ever increasing degree. 

The matter has been investigated theoretically by F. Weiden- 
hammer (litt. 5, 1955) and some hints of it are also found in 
litt. 1-4, but only considered as a question of "stability", not 
realizing that it is a matter of coupled subharmonic vibrations, 
In usual engine plants the phenomenon would only occur at 
impossible high numbers of revolution, but theoretically it is 
very interesting. 

None of the vibration components has an even number of 
vibrations pr. rev. but together they have 2 ~/ rev. At the 
main position the mean value of the phase angles of the two 
components will at resonance be «/4. If the A-vibration at 
this position has a phase angle 9, , the B-vibration will have 


phase angle = ; - Qa. Then the resulting fictive force in 
the A-vibration has phase angle 


leading the B-vibration by 3 » and inversely the fictive force 


{ 
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in the will have phas¢, leading the 
at 
2 


tion by $ - Outside of resonance will approximate 3 


or 0 (resp. above and below resonance), and the action will 
remain completely reciprocal. 

The conditions for steady vibrations at resonance and with 
equilibrium value of the damping will be: 


q.°2 w 
4 a a 1 a 
ae B — = Pp 
b 
-2 

1 qa 4 
which combine into 

4 
Pa Py = (47) 


Usually pa and pp, respectively Bz) and Ppa will be rather 
similar, and (47) is tnen analogous with ce a), sect. 8. Fora 
damping smaller than this the range of increasing vibrations is 
given, analogous with (12a), approximately as 


Litt. 5 has a very similar result found from the stability 
criterion through complicated matrix calculations and confined 
to the special case of a 4-cylinder engine. The method does 
not seem suitable for generalization, and when - as in the 
present statement - a generalization is made, matrix calcu- 
lation becomes obsolete. It is, after all, nothing but vibra- 
tion analysis carried out in a very obscure manner, 

At resonance with equilibrium damping the damping work 
pr. sec. will be the same in both vibration forms, and with 
smaller damping the steady vibrations at the points (48) will 
have the same kinetic energy in the two vibration forms. 

The amplitudes will be approximately inversely proportional 
to the frequencies, 

Now, as in the case of single subharmonic vibrations, 
we must also consider the cylinder impulse K, which, through 
impulse distortion, will have a similar effect as the 2. order 
mass variation. Perhaps it will be more instructive tq con- 
sider the general case of a cylinder impulse Kp (of compa- 
ratively low order number p) causing coupled subharmonic vi- 
brations at a speed 


oto, 


wW = 


° P 


= 


It would here be convenient to construct a "coupled static 


amplitude" 
(K4,4, ) 
pab 
Yim a, =m 4) 


depending on the same type of vector summation as Bap by (43) 
but with a somewhat different denominator. The equilibrium 
condition corresponding to (47) would then be: 


Vea", + P ab (49) 


also analogous with (11) in sect. 3, Similarly we get for the 
critical range 


E ei Paap)” PoP, (50) 


analogous with (12). 

In the cases p = 1-3 where we have both mass variation 
and impulse distortion we can use formulae (47) and (48) when 
in the calculation of Bab and Pha by (43) and (45) B, is replaced 
b 

y pK 
6 
pe 


or we can use formulae (49) and (50) when in the calculation 
of apab we replace Kp by 


Bo mw 2w 
K pc ca K + K 
P P P po 


that means adding a fraction of the oscillating-mass impulse 
pr. cyl. Ky, (see formula (15), sect. 4. It is positive for 
p = 1, negative for p = 2-3), This fraction was in sect. 8 
found = 4 in the case of single subharmonic vibrations, 
corresponding to Wa = 

Perhaps the best way is to define the coefficient of 
coupled mass variation as 


8 z 4,4) 
pab 2 2 
Vz m 4, m 4) 
with the full value of 6,,, and retain as calculated for 
the cylinder impulse only. (49) will then be: 


4 4 


analogous with (11b), sect. 8, and we get a similar extension 


of formula (50). 


It might be expected that there would be possibility of 
coupled subharmonic vibrations also at the speeds 


but it is not so. In these cases we have 


= 4,7 P >0 


therefore ° or 
% = 4, +P 


There is no shift of sign in the q-values and neither in the 
phase angles, consequently no possibility of phase adjustment. 
We have therefore in these cases no excitation of subharmonic 
vibrations, only a "stiff coupling'' working in that way, that 
if there are free vibrations and the damping is small the vi- 
brations will shift periodically from one vibration form to the 
other, with bigger amplitude at the lower frequency. If one 
of the vibration components is excited by an outer impulse 
(which, however, would require q-values as whole numbers) 
the other component is also excited, and the damping coeffi- 
cient will be increased in the proportion 


4 4 2 


Pa PL 


4 


It would not seem improbable to find such cases of dam- 
ping increase in practice, perhaps specially in the case 
wy = 2 when a big impulse of order p would be in direct 
resonance with vibration form A and at the same time the 
artificial impulse of order 2p in resonance with vibration 


form B. The excited vibration B would then work back as 
damping on vibration A. 
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It is now seen, that the vibration analysis for linear multi- 
mass systems given in section 2 is valid also for non-linear systems, 
namely for the system linearized and submitted to fictive forces. 

The only difference is, that in the non-linear system there are 


coupling effects between the various vibration forms, which are not 


present in the linear system without damping. 

It also follows that even in the linear system there may be 
similar coupling effects due to an unequally distributed damping. 
We can namely consider the damping forces as fictive forces on a 
damping-free system, Produced f.inst, by the main vibration in 
vibration form A they are all in the same phase, and their effect 
on vibration form B will depend on the numerical summation 
zra4. But zm4,4, = 0, consequently there will be a coupling 
effect if the damping is not finely distributed with a proportional 


amount on each mass in the system, 


41. Calculations in the actual case 


Due to the symmetry 2408 is 0 in all vibration forms, and 
we have to do only with direct resonance of 8. order and secondary 
resonance of 6. order. 

In vibration form A, fig. 5, we have Pan = 0,11 and 


0.59- 1074 


i 6.25+ 10 k = 0,094 bs 


direct resonance: 0.12° 


The correction member in (41 a) is here found as as, = 0.001, 


negligible in comparison with $8 3 = 0, 041. 


From vibration form B we find By 2" 0.09 and 
* 0.21+-10 0.26+ 10 dg = 0,26 0. 09 
0.01° 
From vibration form C: e. = 0,15 
0.17+- 410 e,= 0,19- 10 e, = 0.19 0.15 


Fig. 10 shows the principle of finding phase angles, when all 
impulse summations =(4) are made in the usual way for crank 1 
in top. The vectors are rotating same way as the crankshaft, and 
the diagrams are then identical for both directions of revolution, 
except that the initial phase angle Y of the cylinder impulse (com- 
mon for all cylinders) should be added in the direction of revolu- 


tion, 

The forces from vibration form B and C are inconsiderable 
in this case, Still the investigation of phase angles shows, that 
all 4 forces are nearly in phase together. The resultant is 


-4 
0.01+ 10 


@« 


Fig. 10. Phase angle for fictive force. 
0.72: 10°* corresponding to a static stress in the shaft at the node 
point T, = 3,80 kg/cm’. With a damping coefficient p = 0,020, 
usual for engines of this type, we finally get a vibration stress 


maximum = 190 kg 


in close agreement with the measurements both at and about the 


resonance point. 


Crank position -diagram K,- diagram 
1 (4) 
1 / 1 1 
j 10 bs 
/ 
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42. Conclusions 


The common, linear theory of torsional vibrations, based on 
differential equations with constant coefficients, is supplemented with 
the analysis of the total vibration in its components of the different 
free vibration forms of the system, which analysis is indispensable 
in any higher vibration theory, especially in cases concerning non- 
linear vibrations in multi-mass systems, As this analysis does 
not seem to be known by the authors of handbooks on vibration 
theory, it was found necessary here first to give a short descrip- 
tion of it and its terminology. 

The ''second order", non-linear theory of crankshaft vibra- 
tions, based on differential equations with variable coefficients, com- 
prises the effect of mass variation due to the oscillating masses, 
and of impulse distortion due to the vibrations. 

Two interesting phenomena are caused by these effects which 
could not be explained by the simple, linear theory, namely sub- 
harmonic vibrations excited by impulses or mass variation compo- 
nents of frequency double of that of the vibration, and "secondary 
resonance" by which resonance vibrations are excited by impulses 
in resonance with a secondary component of the free vibration, of 
frequency different from the main vibration. Such secondary com- 
ponents (in the same vibration form) are caused by the mass va- 
riation. 

These vibration phenomena are here calculated by means of 
“artificial impulses" and "fictive forces" representing the two 
above mentioned effects, whereby the crankshaft system can be 
regarded as a linear system on which the fictive forces act as 
additional outer impulses, The method is a "rationalization" of the 
mathematical method of iteration and gives a great simplification 
of detail calculations, especially in multi-mass systems, There is 
a close connection between the before mentioned vibration analysis 
and the application of fictive forces, because these forces are pro- 
duced by, and acting upon, the "pure'' vibration components in the 


the 


various free vibration forms, and each of these can be calculated 
exactly like one-mass systems. It is only by combination of the 
two matters (analysis and fictive forces) that they come to their 
full advantage, and this is perhaps the reason why none of them 
has so far been admitted to the handbooks, treating one-mass 
problems with great subtlety but cautiously refraining from any 
serious attempts on multi-mass systems, 

The unexpected vibrations in 10-cylinder engines, mentioned 
in the introduction, were found to be a case of ''secondary re- 
sonance", and the case is fully described and calculated, thereby 
also showing the smallness of certain corrections which may be 
disregarded in practical calculations, 

Secondary resonance is a phenomenon common to all kinds 
of non-linear systems, and the calculation method here advanced 
may be useful in other cases, The principle is, that the impulse 
in secondary resonance will in the first place produce ordinary 
forced vibrations in the system, and these forced vibrations will 
then, through the special mechanism of non-linearity in the system 
considered, produce the exciting force for the main vibration, 
Secondary resonance in one vibration form may also be provoked 
by forced vibrations in other vibration forms of the system. The 
effect is fully investigated but found very small in the present case. 

Subharmonic vibrations are also peculiar to all kinds of non- 
linear systems and due to similar effects as the secondary re- 


sonance - therefore the investigation of the two phenomena cannot 


be separated - but there is in our case a great difference between 
them. Secondary resonance in crankshaft systems is working 

exactly like ordinary resonance, only caused by an impulse which 
would not be expected in the linear theory, but the subharmonic 
vibrations are of a very different nature, arising sharply and 
unlimited by damping, if the damping is below a certain treshold 
value. The exciting fictive force is here produced by the vibration 
itself, not by a constant forced vibration, and the exciting force , 
is therefore increasing with the vibration. In the curious case of 
"coupled subharmonic vibrations" free vibrations in two different 


vibration forms will each produce fictive forces exciting the other 


vibration form, 


j 
3 
as 
: 


In crankshaft systems there will normally be damping enough 
to prevent subharmonic vibrations, or it will in any case be easy 
to provide for it, but secondary resonance must be taken as it is 
and included in the surveyance of all the other factors determining 


crank position and shaft dimensions, 
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